arXiv:astro-ph/9706078vl 9 Jun 1997 


Hydrodynamics of GRB Afterglow 

Re’em Sari 

Racah Institute, Hebrew University, Jerusalem 91904, Israel 

and 

Institute for Advanced Study, Princeton, NJ 08540 


Received_; accepted 




ABSTRACT 


The detection of delayed emission at X-ray optical and radio wavelengths, 
“after-glow”, suggests that the relativistic shell which emitted the initial GRB 
due to internal shocks decelerates on encountering an external medium, giving 
rise to the after-glow. We explore the interaction of a relativistic shell with 
a uniform inter stellar medium (ISM), up to the non relativistic stage. We 
demonstrate the importance of several effects that were previously ignored, and 
must be included in a detailed radiation analysis. At a very early stage (few 
seconds), the observed bolometric luminosity increases as t 2 . On longer time 
scales (more than ~ 10sec), the luminosity drops as t _1 . If the main burst is 
long enough, an intermediate stage of constant luminosity will form. In this 
case, the after-glow overlaps the main burst, otherwise there is a time separation 
between the two. On the long time scale, the flow decelerate in a self similar 
way, reaching non relativistic velocities after ~ 30 days. Explicit expressions 
for the radial profiles of this self similar deceleration are given. Due to the 
deceleration and the accumulation of ISM material, the relation between the 
observed time, the shock radius, and its Lorentz factor, is given by t = R/l&^c 
which is a factor of 8 different from the usual expression. The majority of 
particles are those of the original ejecta (and not the ISM) up to about ~ 900s. 
These particles reach sub-relativistic velocities on a time scale of ~ 2hours, 
well before the flow becomes sub-relativistic. Therefore the ejecta particles are 
probably unimportant for most of the after-glow radiation. We show that even 
though only a small fraction of the energy is given to the electrons, most of 
the energy can be radiated over time. If this fraction is greater than ~ 10% 
radiation losses will significantly influence the hydro dynamical evolution at early 
times (less than lday). 
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1. Introduction 

The isotropy of GRBs angular distribution combined with the non-homogeneous 
distribution suggests that GRBs originate from cosmological distances and therefor radiate 
energies of order of 10 51 erg. Considerations of optical depth then show that the bursts 
are produced by the dissipation of kinetic energy of highly relativistic shells with Lorentz 
factor r) > 100 (see Piran 1996 for review). This dissipation can be either due to internal 
shocks or due to the surrounding ISM. Sari & Piran (1997) have shown that deceleration on 
the ISM could not give rise to the variability observed in the bursts (unless the process is 
very inefficient and involves much more than 10 51 erg), while internal shocks could produce 
efficiently the observed fluctuations (Kobayashi, Piran and Sari 1997), if the “inner engine” 
has considerable fluctuations. It is therefore likely that the main GRB is due to internal 
shocks. 

These cosmological models predict that after the main GRB event the ejecta decelerate 
due to interaction with the ISM, emitting radiation at longer and longer wavelengths 
(Paczynski & Rhoads 1993, Katz 1994, Meszaros & Rees 1997). This emission has been 
detected recently for several GRBs due to an accurate determination of their position. The 
quantitative agreement between the deceleration models and the measurements is good 
(Waxman 1997a, Wijers, Rees & Meszaros 1997, Waxman 1997b). Since the quality of data 
for the after-glows is higher than for the burst itself, more quantitative results are needed. 

In this letter we explore the hydrodynamics of the deceleration of a relativistic fireball 
on a uniform ISM. We discuss the relation between the rise time of the after-glow and the 
time of the main burst. We use the analytic solution found by Blandford and McKee (1976) 
to describe the swept up ISM and derive an expression for the position of the original ejecta 
and its Lorentz factor. We show that although only a small fraction of the internal energy 
is given to the radiating electrons, a considerable amount of the energy can be radiated over 



the deceleration period. 


2. The Rise of The After-Glow 

The problem of deceleration of a relativistic shell onto the ISM is determined by four 
parameters: the initial shell’s Lorentz factor rj, the energy of the shell E = i?52l0 52 erg, the 
width of the shell (observer frame) A and the ISM density n = n jcrri -3 . The basic details of 
the interaction between the shell and ISM where given in Sari & Piran 1995 and we briefly 
review the main ideas. The treatment in this section is approximate and correction factors 
of order unity may need to be included in a more precise treatment. 

When the shell encounters the ISM, two shocks are formed: a forward shock accelerating 
the ISM and a reverse shock decelerating the shell. The forward shock is always highly 
relativistic since the initial Lorentz factor rj 1. Let / be the density ratio between the 
shell and the ISM given by 

/ = E/I6nrj 2 Anm p c 4r yH 2 . (1) 

where we used R = 2cy 2 t for the radius of the shell, t is the observed time. The reverse 
shock is relativistic if / < rf reducing the Lorentz factor of the shell to 7 = ^V 2 / 1 / 4 /-^2 
and Newtonian if / > if making only negligible change to the shell’s Lorentz factor i.e., 

7 = 77 . At early stages, f » rj 2 so the reverse shock is Newtonian and the shell’s Lorentz 
factor equals its initial value 7 = 77 . However due to the increase in the area of the shell, it 
produces internal energy in an increasing rate of 

L = 327rc 5 mn P 7 8 t 2 = 2.5 x 10 5 ° 7 f 00 7 ii £ 2 ——, (2) 

s 

where we use t s for the time in seconds. Assuming that the cooling is fast (Sari, Narayan 
& Piran 1996), the observed bolometric luminosity is proportional to the internal energy 
increase rate and is therefore also given by Eq. (||). This behavior will continue until either 
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the shell has given the ISM an energy comparable to its initial energy at 
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whichever comes first. As in Sari & Piran (1995) we define the ratio between the two 
expressions as 
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For £ > 1, the energy is dissipated to internal energy before the Lorentz factor of the 
shell is reduced considerably. If £ < 1, then the reverse shock turns relativistic before the 
kinetic energy of the shell was emitted. In this case the Lorentz factor decreases with time 
according to 

7(0 = y /2 /(i) 1/4 /V2 = 300 Eg ( 6 x 101 , cm ) n^h;V‘. (6) 


Note that at this stage the Lorentz factor 7 is independent of its initial value 1 /. Substituting 
this in Eq. (® we get that the luminosity is constant in time and given by: 


L = 


E 

2 A /c 


(7) 


This stage will continue until the shell has given the shocked ISM energy comparable with 
its own at 

t E = 2A/c. ( 8 ) 


At this time the shell has decelerated to Lorentz factor of 
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(9) 
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independent of the initial Lorentz factor 77 . After the time given by Eq. (|3|) or ^ the 
flow will be described by a self similar solution as the ISM energy is now constant and 
comparable to the initial energy of the shell E. As we shall show in the next section, from 
this time 7 oc R ~ 3 / 2 oc t ~ 3 / 8 and therefore the observed luminosity decreases as £ _1 . This 
behavior is illustrated in figure 1 . 


3. Relation with the main burst 

If the main burst is produced by internal shocks, then the width of the shell, A, can 
be inferred directly from the observed main burst duration A = ct m b . For long bursts, 
t m b ~ 20 sec and A ~ 6 x 10 11 cm while for short bursts t m b ~ 0.1 sec and A ~ 3 x 10 9 cm. 
The initial Lorentz factor must satisfy 77 > 100 for the emission of the main burst not to be 
opaque. The reverse shock is therefore likely to be Newtonian for short bursts and might 
be relativistic for long bursts. Both cases are therefore of physical interest. 

If the reverse shock is relativistic, then the observed peak of the after-glow emission is 
flat and overlaps the observed GRB emitted by internal shocks. Both end after an observed 
time of ~ A/c. If the reverse shock is Newtonian, then the after-glow peaks on given by 
Eq. (3j) which is longer by a factor of 

> 1 ( 10 ) 

than the main burst duration A/c. The duration and luminosity of the main burst and the 
after-glow rise are shown for the Newtonian and relativistic cases in figure 1. 

In both cases, the properties of the main burst and the after-glow are very different. 
The main burst is usually highly variable (depending on the internal structure of the shell) 
while the after-glow, which is due to external shocks is expected to be smooth (Sari & Piran 
1997). The after-glow’s spectrum should peak, in the beginning, around 30KeV-10MeV 



t(s) t(s) 


Fig. 1.— The luminosity from the ISM as function of time in the Relativistic (left frame) and 
Newtonian (right frame) cases is drawn in dashed line. At the early stage, the Lorentz factor 
is constant and the luminosity increases due to the increase in shell area. When the ISM 
has energy comparable with the total energy (t — £g), a self similar solution is established 
and the luminosity drops as f -1 . If the shell is thick (typical for long main bursts) the 
reverse shock becomes relativistic at tjy, before the self-similar solution is established and 
some deceleration begins, leading to constant luminosity. Solid line gives the luminosity of 
the main GRB. Both frames use E = 10 52 erg and 77 = 300. The behavior L ~ t _1 and 
7 ~ £ -3 / 8 continues up to the non-relativistic stage which is about 30days. 
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depending on the fraction of internal energy in electrons and magnetic field (Sari, Narayan 
and Piran 1996). IF the peak energy is too high it might not be observed in the first 
stage by the BATSE equipment. However later as the ejecta decelerates the emission peak 
decreases in time and should cross the soft 7 -ray region. 


4. Self-Similar Solution and Properties of the Original Ejecta 


We begin with a simple consideration based upon conservation of energy. When most 
of the energy has been given to the ISM, and assuming that radiation losses are small, the 
energy in the shocked ISM is constant and approximately equal to the initial kinetic energy 
E. The shocked ISM rest mass is M oc R 3 . Since it was heated by a relativistic shock its 
energy in the observer frame is ~ M 7 2 . Comparing this with the constant total energy of 
the system E we get that 

7 oc iT 3/2 . ( 11 ) 

This is also the scaling law for the shock wave Lorentz factor T since for a relativistic shock 

r = ^7- 


The scaling law, Eq. ©, implies a quantitative but important change in the relation 
between t, R and 7 . Photons that were emitted from the shock while it has propagated 
a small distance SR will be observed on time-scale of 5t ~ 5R/2Y 2 c. Integrating this over 


time using the scaling law |TT] we get t = R/8T 2 c, or 

R 
t = 


167 2 c 


( 12 ) 


Compared with the commonly used expression i?/2y 2 c (Meszaros and Rees 1997, Waxman 
1997a, Wijers, Rees and Meszaros 1997, Waxman 1997b), this expression is factor of 8 
smaller. This difference is important when trying to fit quantitatively the observed after¬ 
glow data. Note however that the differential relation is independent of the deceleration, 
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and is therefore given by St = SR/Ary 2 c. 


Blandford and McKee (1976) have described an analytical solution for the case in 


which the scaling law [XT] applies. Using their solution with several simplifications and some 
algebraic manipulations we get 


n(r , t) 
7 (r,t) 
e(r,t) 


= 4nj 1 + 167 2 (1 — r/R) 


-5/4 


= 7 


1 + 167 2 (1 — r/R) 


- 1/2 


= 4 nm p c 2 7 2 1 + 167 2 (1 — r/R) 


1 — 17/12 


(13) 


where n(r,t), 7 (r,t) and e(r,t) are, respectively, the density, Lorentz factor and energy 
density of the material behind the shock (not to be confused with the ISM density n and 
the Lorentz factor of material just behind the shock 7 (t) = 7 (R,t)). The scaling laws of 
R(t) and 7 (t) can be found using these profiles and demanding that the total energy in the 
flow would be equal to E: 

17Et \ 1/4 


m = 


^ 7 Tm p nc 


J 


= 3.2 x 10 16 E^ 4 n/ 1/4 ty 4 cm 


**> = y^r= 26oE * ws 


(14) 


This solution can serve as a starting point for detailed radiation emission calculations and 


comparison with observations. The scalings [14] are, of course, consistent with the scalings 
0 and [T| which were derived using conservation of energy, but supply the exact coefficient 
that can not be produced otherwise. The time for which the flow behind the shock becomes 
sub-relativistic follows from Eq. fll4]) as 


t = 30E , 52 3 ti 1 1//3 days. 


(15) 


We now turn to see the role of the particles from the original ejecta in the flow. The 
number of protons (or electrons) in the ejecta is E/rjm p c 2 , and is larger than the number of 
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ISM protons swept by the shock ( 47 ri? 3 n/ 3 ) up to the time 

/ q4 E \ 


t = 


= 850 E 5 A 1 1/3 )?3to /3 s 


(16) 


y 17 3 4 4 7 t nm p r) 4 c 5 j 

We now use the solution [13] to determine the evolution of the original ejecta material (or 
any other fluid element). Taking the derivative of the fluid Lorentz factor along its line of 
motion we get 


dru _ _l7e 

dt 8 t 


(17) 


so that a fluid element which had a Lorentz factor 70 at to 5 will have 

/ t \ “ 7 / 8 

7e ( t)=7°y • (18) 

The exponent —7/8 shows a fast deceleration relative to the shock deceleration exponent 
—3/8. For the Newton case, the ejecta had a Lorentz factor 7 at the beginning of the 
self-similar stage at time T E given by Eq. ([5|), therefore reaching sub-relativistic velocities 
after 

tar = ^h 8/ ' = 0-9 E 52 3 ni 1/3 %00 /21 llours ( 19 ) 


For the relativistic case we get 

/ A \ 4 / 7 

tsr = t E l* / 7 (t E ) = 2.6 ElJ 2 7 r h l/7 ( 6 x 10 n cm j hours - ( 20 ) 

This time scale of ~ 2hours is much shorter than the time in which the whole flow becomes 
relativistic ~ 30days. It is therefor expected that most of the afterglow radiation, that was 
seen on time scale of days and even month, is not related to the particles of the initial 
ejecta but to the shocked ISM. 


5. Radiative Corrections 

I 11 the previous sections we have assumed that the energy in the system is constant. 
This assumption can not be strictly correct since the radiation takes some energy from the 
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relativistic shell. We define e e to be the fraction of the internal energy that is radiated, and 
lost from the system. Typically this should be the fraction of energy given by the shock to 
electrons and is estimated to be around 10% (Waxman 1997a,b). This number seems to 
be negligible, and therefore the energy loss was neglected by previous analyses. However 
the deceleration occurs over several orders of magnitude in time and Lorentz factor and the 
fireball energy E is given again and again to newly heated material, leading to more and 
more energy losses. 


The energy loss rate during the deceleration is given by 40967 rc 5 nm p 7 8 t 2 (the coefficient 
is different from Eq. 2, due to the relation R = I 67 2 ct which is relevant in the deceleration 
stage) multiplied by e e . Substituting the expression for y(t) from the self similar solution 
(Eq. [L4|) we get 


dE _ 17 E 

dt 16 e t ’ 


( 21 ) 


so that 


/ + \ - 17 V16 

E(t) = E 0 (-) • (22) 

Since the observed initial time of the afterglow is about 10s then after about a week the 
energy is reduced by a factor of ~ 3 for e e = 0.1 or a factor of ~ 30 if e e = 0.3. These 
factors must be taken into account given the accuracy of current data. 


The derivation of the above exponent, —17e e /16, used the exact coefficients in Eq. 
(O). which were obtained from the self-similar solution. Without this solution the exponent 
could only be estimated approximately. Note that the use of Eq. ( |T4| ) is valid as long as e e is 
small enough that the energy loss could be approximated as a small “radiative correction”. 


This radiation losses will also slightly effect the scaling of the shock radius and Lorentz 


factor as function of time. The approximate scaling including the radiation losses can be 
obtained by substituting Eq. (|2^) into Eq. ([14]). 
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6. Discussion 

We have explored the early evolution of the interaction of a relativistic shell with the 
ISM. If the main GRB is short enough, separation is expected between the main burst 
and the afterglow luminosity peak, while if it is long enough an overlap is expected. This 
property might be detectable in BATSE’s data. 

We have used the self similar solution derived by Blandford and McKee (1976) to 
obtain an explicit expression for the radial profile in the self similar stage. This solution can 
be used in further analyses when considering a more detailed calculation of the radiation 
from the heated ISM. 

A relation between the shock’s radius, the material Lorentz factor and the observed 
time was found to be t = R/1 Q^ 2 c instead of the commonly used expression t = R/2^ 2 c 
due to the fact that ISM is collected so the shock moves faster than the material behind it, 
and due to the deceleration of the shell, having higher Lorentz factor at earlier time. This 
relation was obtained assuming that the radiation is emitted from the shock front. On long 
time scales ~ lday, when the cooling of electrons is not fast enough it might be that the 
width of the radiating zone will smear the observed radiation over longer time scales. 

The role of energy loss due to the radiation was found to be non-negligible even if the 
part of the internal energy that is radiated at each time is small. Thus radiation can reduce 
the total energy in the system after a week by a factor of 3 if e e = 0.1 or a factor of 30 if 
e e = 0.3. 

The data of GRB970228 and GRB970508 fit radiation models, with in a factor of 2, 
without taking into account radiation losses (Waxman 1997a,b). We can therefor roughly 
estimate e e < 0.1. On time scale which is more than ~ lday, the electron’s cooling time 
is long so only a small fraction of their energy is radiated. Since most of the observations 



where made after ~ lday, the fraction of the energy that is given to the electrons can be 
high (more than 10%) without leading to considerable energy loss, and with no effect on 
the observations made after ~ lday. ffowever, in such a case, the energy losses at earlier 
time will be considerable and will therefor require a much higher initial energy. 

The author thanks The Institute for Advanced Studies for warm hospitality and Eli 
Waxman, Pawan Kumar, John Bahcall, Tsvi Piran, Jonathan Katz, and Shiho Kobayashi 
for helpful discussions. 
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